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Abstract 

We have simulated the response of a high energy neutrino telescope to the 
stream of low energy neutrinos produced by a supernova. The nominal threshold 
of such detectors is in the GeV energy range. The passage of a large flux of 
MeV neutrinos during a period of seconds will nevertheless be detected as an 
excess of single counting rates in all individual optical modules. Detectors under 
construction, which consist of roughly 200 modules, will be able to detect a 
galactic supernova at or above the 5 a level. The rate of fake signals is, however, 
too large for the telescope to serve as a neutrino watch. Such capability requires 
detectors with roughly 3 times the number of optical modules, thus within easy 
reach of the next generation detectors. 



The neutrino events detected in the Kamiokande[|l|, IMB|0], Baksan[Q and LSD[0 de- 
tectors prior to the optical discovery of supernova 1987A represented a most remarkable 
birth of neutrino astronomy. Despite observation by four experiments, the data has left 
us with a variety of questions and uncertainties. Most prominent is our inability to un- 
derstand the time of the Mont Blanc neutrino burst [^] and the directionality of the 1MB 
events [^]. This underscores the importance of collecting as much information as possible 
when presented with the rare opportunity of observing the next nearby supernova. We here 
show that high energy neutrino detectors, presently under construction 0, can observe the 
neutrino bursts from galactic supernovae even though they have nominal thresholds of GeV 
energy, i.e. almost three orders of magnitude above the MeV energy of supernova neutrinos. 
Simultaneous detection of supernova neutrinos, together with other devoted experiments|P], 
will help overall statistics. Good relative timing between different experiments may allow 
the determination of the location of the supernova if it is shielded from us by large amounts 
of matter. 

First generation high energy neutrino telescopes consist of approximately 200 optical 
modules (OM) deployed in deep, clear water or ice shielded from cosmic rays. Coincident 
signals between the OMs detect the Cerenkov light of muons with energy in excess of a 
few GeV. Also electromagnetic showers initiated by very high energy electron neutrinos are 
efficiently detected. The idea has been debated for some time whether these instruments 
have the capability to detect the neutrinos from a supernova despite the fact that they are 
in the MeV range. The production of copious numbers of positrons of tens of MeV energy in 
the interaction of z/g with hydrogen, will suddenly yield signals in all OMs for the 10 seconds 
duration of the burst. Clearly such a signal, no matter how weak, will become statistically 
significant for a sufficient number of OMs. We here perform a complete simulation of the 
signal and its detection and conclude that the 200 OMs of detectors such as DUMAND and 
AMANDA are sufficient to establish the occurence of a neutrino burst in coincidence with 
the optical display of a supernova. We also show that the same detectors can actually serve 
as a supernova watch, i.e. a "fake" signal occurs less than once a century, by increasing the 
number of OMs by a factor of three. This is much less than the roughly 7000 OMs which 
are projected for a next-generation detector 0. 

Although aspects of the observations of SN1987A left some lingering doubts about su- 
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pernova models|^, they remarkably confirmed the estabhshed ideas for the supernova 
mechanisms 0. At collapse the core is expected to release energy in a prompt z/g burst 
lasting a few milliseconds. Most of the energy is however liberated after deleptonization 
in a burst lasting about ten seconds. Roughly equal energies are carried by each neutrino 
species. The time scale corresponds to the thermalization of the neutrinosphere and its 
diffusion within the dense core[^. Since the v^. cross-section [^] for the inverse beta decay 
reaction on protons exceeds the characteristic cross sections for the other neutrino flavors, 
Vf. events dominate by a large factor after including detection efficiency. In this reaction 
free protons absorb the antineutrino to produce a neutron and a positron which is isotrop- 
ically emitted with an energy close to that of the initial neutrino. For the purpose of this 
letter we use typical parameters, derived from SN1987 observations, which are consistent 
with those previously estimated in supernova models. From the energy distributions of the 
observed events the average temperature of the neutrino sphere in SN1987A was deduced to 
be 4.0 MeVH. 

Before discussing our detailed Monte Carlo simulation we present a back-of-the-envelope 
derivation of our final result. After convoluting the 4 MeV thermal Fermi distribution of the 
neutrinos with a detection cross section rising with the square of the neutrino energy, one 
obtains an event distribution peaked in the vicinity of 20 MeV. The tracklength of a 20 MeV 
positron is roughly 10 centimeters and therefore over 3000 Cerenkov photons are produced. 
This number combined with a typical quantum efficiency of 25 % leaves 800 detected photons 
in each event. Assuming these are emitted within vr steradians, i.e. roughly the size of the 
Cerenkov cone, the detection probability becomes a simple function of both the module 
collection area Am and the distance to the positron shower R: 
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For < R^ (~ 250 Am) the phototube will most likely trigger on the positron while for 
larger R the probability diminishes rapidly. By evaluating the effective volume within a 
cone of vr steradians with vertex at each OM, we obtain an approximate expression for the 
detection volume associated with OMs: 

1 rRat 



Veff ~ P{R)R'dR 
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Here we have integrated up to the attenuation length Ratt of the medium and divided by a 
factor of 4 in order to average over all orientations of the Cerenkov cone. We conservatively 
assumed that the OM has 27r acceptance. It is interesting to note that the effective volume 
is proportional to both the collection area of the OM and the attenuation length. It is quite 
insensitive to the solid angle over which the Cerenkov photons are distributed (vr). For OMs 
such as those used in the AMANDA detector with collecting area Am = 0.028 m^ and an 
attenuation length Ratt = 25 m typical of ice[|lT|, we obtain Veff ~ 130 m^. This result can 



be used to rescale SN1987 observations to a supernova at a distance d^pc- From 11 events 
observed in 2.14 kton Kamiokande detector we predict: 
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for a detector with Nm optical modules. For a 130 m'^ effective volume of each of the 200 
OMs we obtain 5300 events. 

We now require a meaningful detection of this signal in the presence of the continuous 
background counting rate of all phototubes. Over the 10 s duration of the delayed neutrino 
burst from a supernova, the rms fluctuations of the combined noise from all the OMs is: 



CTip.e. = Z/ip.e. Nm (4) 

where the background counting rate in each module at the 1 photoelectron level is represented 
by z^ip.e.- The probability that the noise in the OMs fakes a supernova signal can be estimated 
assuming Poisson statistics. The expected rate of supernova explosions in our galaxy is 
about 2 X 10^^ y^^. If the detector is to perform a supernova watch we must require that the 
frequency of fake signals is well below this rate. The signal should therefore exceed n„ > Q 
which corresponds to a probability of 9.9 x 10~^°. The corresponding number of 10 seconds 
intervals indeed exceeds a century. Clearly the requirement can be relaxed if we just demand 
that the detector can make a measurement in the presence of independent confirmation. For 
an average noise rate of 1 kHz, a typical value for the OMs in AMANDA, the rms fluctuation 
of the 20 million hits expected in an interval of ten seconds is 1400. This implies that 



detection of a galactic supernova is near the 4 a level for the 200 module configuration, while 
detection should not represent a problem for the next generation detector which consists of 
7,000 OMs. Since the signal in the present detector is marginal, it is necessary to do a more 
realistic calculation of the event rate. We will conclude that our rough estimate is somewhat 
conservative. 

Background noise in the modules clearly plays a critical role so that low noise environ- 
ments such as ice have an intrinsic advantage. Furthermore, the noise is expected to be 
reduced drastically at the 2 photoelectron level. This will unfortunately also imply a re- 
duction in effective volume for event detection as we will see further on. Signal to noise is 
proportional to the ratio K/z/'sA'- Obviously increased attenuation length in the medium 
and larger effective area of the OM results in an enhanced effective volume. Considering 
parameters appropriate for DUMAND, an attenuation length of 40 m in water and OMs 
with double diameter, we expect a factor 10 increase in effective volume per optical module. 
This should readily compensate for a noise rate higher by a factor 100. We therefore expect 
DUMAND and AMANDA to have comparable sensitivity as supernova detectors. 

For a complete calculation we have combined a detailed electromagnetic shower Monte 
Carlo, initially developed to evaluate the radio emission by cascades in ice[|^, with the 
AMANDA Monte Carlo. The shower Monte Carlo is a fast three dimensional routine which 
simulates the dominant low energy processes: M0ller, Bhabha and Compton scattering as 
well as electron-positron annihilation, continuous energy loss and multiple elastic scattering 
as well as the bremsstrahlung and pair production processes which dominate at high energy. 
We added the capability to simulate the emission of Cerenkov photons by cascade particles. 
Our event file consists of over 10,000 events sampled from a 4 MeV temperature Fermi-Dirac 
neutrino distribution weighed by a cross section rising with the square of the neutrino energy. 
Photon detection is simulated using the AMANDA Monte Carlo to correctly account for the 
effects of attenuation in deep polar ice, optical module efficiency as a function of photon 



wavelength and detector geometry] 13 



The final result can be quoted as an effective volume for the entire detector of 23,000 m'^i 
2, 000 for the first stage AMANDA configuration of 10 strings arranged in a nine 30 m side 



polygon with one at the center. Each string has 20 modules spaced at 10 m intervals [IJ]. We 
therefore obtain an effective volume of 115 m^ per module, close to our crude estimate. Single 
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OMs have a lower threshold than the Kamiokande and 1MB experiments which reconstruct 
Cerenkov cones. This is taken into account by correction factors multiplying the event rates 
to be entered in Eq. (3), which we evaluated to be 1.4 (5.4) for Kamiokande (1MB). We 
thus obtain 7700 (9070) events in the 200 modules of the 9+1 configuration of AMANDA. A 
trigger could be implemented by monitoring the sum of all singles rates in ~ 1 s. intervals, 
and requiring positive fluctuations in this sum for several such intervals. 

At the 2 p.e. level we obtain an effective volume of 5,900 m^ ± 1, 400. If the OM's noise 
rate is reduced by over a factor 16, signal-to-noise will actually be improved by working at 
the 2 p.e. level. Clearly the statistical significance of the signal will be improved by including 
information from higher level triggers. 

Typical noise rate for AMANDA modules imply that aip_e. ~ 1400. The calculated event 
rates for a supernova bursts in the center of the galaxy therefore correspond to a 5.4 (6.4) 
sigma effect when rescaling the observed signals from SN1987A at Kamiokande (1MB). This 
should provide a sufficiently clean signal. We can combine Eqs. (2) and (3), including the 
correction factor, to obtain the expected signal: 



n„ = 0.35 X N, 



102 ' ^ - , r. , n2 



125 m3 



8 kpc 

dkpc 



(5) 
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We here scaled to the Kamiokande event number. The constant is 0.42 for 1MB. The re- 
sults are encouraging since a AMANDA supernova watch can be performed at the 6-sigma 
level with only 340 (240) optical modules according to Kamiokande (1MB) observations of 
SN1987A. A next-generation detector with over 7,000 modules should provide a sharp signal 
for a supernova at the galactic center and should operate as a supernova watch to twice the 
distance to the galactic center, i.e. covering all the galactic disk. 

We end with a warning about the statistics. Clearly in a realistic analysis penalty factors 
will be associated with various trials made to identify a burst. This will, however, not alter 
our positive conclusions. A modest increase the number of OMs can absorb the effect of a 
large number of trials, e.g. associated with a sliding window to identify the 10 seconds burst. 
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